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Polarized proton collisions will be studied at RHIC up to a total center of mass 
energy of 500 GeV, starting in 2002. An overview of the RHIC-spin program, and 
the critical components of the PHENIX and STAR detectors for spin experiments, 
is presented. The premier experiment within the RHIC-spin program is the deter- 
mination of the fraction of the proton's spin carried by gluons. A detailed analysis 
of how accurately this fraction can be determined by the STAR experiment at 
RHIC is presented. 



1 Introduction 

The purpose of this workshop is to estabhsh the primary physics goals of a po- 
larized electron-polarized ion colhder (EPIC) that would become operational, 
at the earliest, in 2005. To predict what the most interesting questions will be 
in the future, it is necessary to understand what questions are being addressed 
in the present, and to establish the expected quality of the answers. With that 
perspective, I will describe plans for an experimental program at Brookhaven 
National Laboratory that will study the high-energy collisions of polarized pro- 
tons, accelerated in the Relativistic Heavy Ion Collider (RHIC). There is great 
expectation within the community that the RHIC-spin program will provide us 
with important information about the spin structure of the proton. Arguably, 
the primary goal of the RHIC-spin experiments is to determine the fraction 
of the proton's spin carried by gluons (the integral AG, defined below). Es- 
tablishing the degree of polarization of the glue within the proton is the next 
essential step in understanding the spin structure of the proton. Under the 
assumption that experiments establishing the gluon polarization are the most 
important, most of this talk will be devoted to that subject. 



2 An overview of the RHIC spin program 

Before delving into the details about measuring the polarization of the proton's 
glue, it is important to give a broader overview of the physics of RHIC spin. 
To date, the highest energy achieved in a synchrotron for a polarized proton 
beam is 24.6 GeV in the AGSld. The RHIC rings will provide polarized proton 
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beams up to a maximum energy of ^250 GeV, an order of magnitude higher 
in energy. Even more impressive is that the RHIC-spin program will provide 
total center of mass energy {^/s) for collisions between polarized protons up to 
500 GeV. Furthermore, the expected luminosity for p + p collisions at RHIC 
(2 X 10^^cm~^^s~^ at ^/s = 500 GeV), will enable a systematic study of large 
transverse momentum (px) processes, where perturbative QCD has been suc- 
cessfully applied to explain much of the data obtained at unpolarized pp and pp 
colliders. The addition of polarization to high-pT collisions between protons 
at very high energies will provide an important test of QCD. 

With the addition of polarization to colliding proton beams, spin asym- 
metries (A): 

are new observables that can be measured. In Eqn. 1, the beam polariza- 
tion(s) are represented by P and are expected to be 70% for the RHIC-spin 
program. Asymmetry measurements are made by determining how the yield 
for some process varies with the polarization state {N±) of the beam(s) and by 
measuring the polarization of the beam(s). Eqn. 1 assumes equal integrated 
luminosities for the two polarization states of the beam(s). For the RHIC- 
spin program, it is expected that transverse single-spin asymmetries (At) will 
approach zero for sufficiently large pt', non-zero values for At are expected 
only from higher- twist contributions. Parity- violating longitudinal single-spin 
asymmetries (Al) are^xpectcd to be quite large in the production of real 
weak bosons {W^, Z*')aor at sufficiently high px, where virtual weak bosons 
contribute significantly El to the force between the interacting partons. 

Double-spin asymmetries involve differences of yields for a process when 
the initial state protons have the same and opposite direction polarizations. 
The (parity-allowed) longitudinal asymmetries (All) are expected to teach 
us about helicity asymmetry structure functions, similar to those probed in 
polarized deep inelastic scattering. Transverse double-spin asvmmetries (Att) 
may enable measurement of the transversity distributions related to the 
transverse polarization of the nucleon's constituents. 

If QCD survives the onslaught of these new high-energy and high-p^ po- 
larization observables, then polarized proton collisions can be used as a tool 
to better understand the spin structure of the proton. As described in de- 
tail below, All is very sensitive to the gluon polarization for several different 
processes. In addition, parity-violating Al measurements for the production 
of and Z° can be directly related to the polarization of valence and sea 
quarks within the proton 13. Since the valence quark polarizations are well 
determined in polarized deep inelastic scattering, direct tests of the Standard 
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Model are possible. By changing the kinematical conditions, p + p ^ W + X 
can be used to selectively probe antiquark polarizations, thereby providing new 
information about the proton's spin structure. Polarization observables for the 
highest pt processes at RHIC may also provide limits on either the existence of 
new vector bosons or possible quark substructures^, competitive with planned 
measurements at the Tevatron. 



3 Preparing RHIC for polarized proton studies 



The RHIC spin program will present 
the opportunity to study observ- 
ables in polarized pp collisions at 
y/s values that are more than a 
factor ofj-£5 larger than previously 
available a. This frontier of high- 
energy spin physics is made techni- 
cally feasible by significant advances 
in the handling of polarized beams 
in synchrotrons and storage rings. 
The new technique employs so-called 
'Siberian snakes', a Xjoncept intro- 
duced in Novosibirsk □ and first ex- 
perimentalb: verified at the lUCF 
Cooler ringQ. The Siberian snake is a 
means of overcoming the large num- 
ber of depolarizing resonances en- 
countered when attempting to accel- 
erate polarized proton (or other light 
ion) beams to high energies in a cir- 
cular accelerator. For RHIC, heli- 
cal dipole magnets, funded by the 
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Figure 1: Schematic layout of the RHIC ac- 
celerator complex, indicating the critical ele- 
ments needed for the acceleration of polarized 
ions. 

RIKEN institute in Japan, will be used to reverse the direction of unwanted 
transverse polarization components during alternate turns around each ring. 
In the absence of such spin manipulations, it is possible to destroy the beam 
polarization by small perturbations that act in a single turn, with their effects 
adding coherently when the ring is operated near a depolarizing resonance. 

The motivation to build RHIC and its two major detectors, PHENIX and 
STAR, originated from the goal of discovering and studying a new state of 
matter, the quark-gluon plasma^ that is thought to have existed at the earliest 
moments after the Big Bang. It is expected that the quark-gluon plasma should 
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be formed in the ultrarelativistic collisions between heavy-ion beams at RHIC. 
Although the major detectors were not designed for the study of p+p collisions 
up to total center of mass energy equal to 500 GeV, many of their subsystems 
have been subsequently adapted to this task. As a result, PHENIX and STAR 
bring complementary strengths to the RHIC-spin program. 

The PHENIX detector has many relatively small acceptance detectors that 
provide fine granularity and precise particle identification. Of greatest rele- 
vance to the spin program are the photon arms that are centered at midra- 
pidity (I77I < 0.35), and the muon arms, spanning 1.0 < jryj < 2.4. The fine 
granularity of the Pb-glass and Pb+scintillating fiber elements of the PHENIX 
electromagnetic calorimeter will provide the capability to separately identify 
single photons and di-photon pairs, produced in the decay of neutral mesons. 
These detectors will be employed for inclusive photon production in the RHIC- 
spin program. As well, they will be used to study high-pT tt'^ production. The 
muon arms provide excellent particle identification, and will be used for the 
study of production, by observing the /i* daughters of the W . In addi- 
tion, the detection of jJL^ iJi~ pairs with the PHENIX muon arms will enable a 
program of measuring polarization observables associated with vector meson 
production and the Drell-Yan process. 

The STAR detector is intended to provide a more global view of a heavy- 
ion collision. The heart of STAR is a 0.5 T solenoidal magnetic field and a 
time projection chamber (TPC), capable of tracking all of the charged par- 
ticles produced in a central Au-Au collision in the nominal pseudorapidity 
interval, |?7| < 2, with full azimuthal coverage. Multiple layers of silicon de- 
tectors around the interaction point will be used for reconstructing primary 
event vertices, and secondary vertices from strange particles. Of greatest rel- 
evance for the RHIC-spin program are the barrel and endcap electromagnetic 
calorimeters (EMC). Construction of_the barrel EMC is underway and should 
be completed by 2002. A proposals for the endcap EMC is awaiting final 
decision concerning funding. If approved, the construction timetable for that 
detector will be comparable to the barrel EMC project. In comparison to 
PHENIX, the STAR EMC has substantially coarser granularity, but provides 
full azimuthal coverage for the pseudorapidity interval, —1 < 77 < 2, with a 
small gap near 77 = -1-1 for TPC and EMC services. The large acceptance 
of STAR makes it ideal for jet reconstruction. As described below, detection 
of photon -t- jet coincidences provides STAR with unique capabilities in de- 
termining AG. In addition, measurements of di-jet production at STAR will 
provide an important cross check of the gluon polarization measured in direct 
photon production. Finally, the STAR EMC will also enable the study of 
production, by observing the daughters of the W . 
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As of July 1999, the status of the RHIC accelerator is as follows. Commis- 
sioning of Au beam in the RHIC rings began this summer. The initial studies 
of Au-Au collisions will commence in November 1999. The first helical dipole 
magnets have been successfully produced, and the initial commissioning of a 
polarized proton beam is scheduled for 2000. Initial studies of p + p collisions 
at low luminosity are planned for 2001. It is expected that full luminosity 
p + p collisions will begin in 2002. An agreement between RIKEN and BNL 
will allow the RHIC-spin program to run for 10 out of the 37 weeks of annual 
RHIC operations. It is projected that at full luminosity, a 10-week run at 
\/s=200(500) GeV will result in an integrated luminosity of 320(800) pb'^ 

4 Methods of determining AG 

After this brief overview of the spin-physics program and the tools to carry it 
out, I'll focus the rest of this writeup on the determination of the fraction of the 
proton's spin carried by gluons, which is expected to be the most important 
result forthcoming from the RHIC-spin program. That fraction is equal to 
twice the integral of the gluon helicity asymmetry distribution {AG{x, Q^)): 

AG(Q2)= / AG{x,Q^)dx= [ [G+{x,Q^)-G-{x,Q^)]dx, (2) 
Jo Jo 

and is a function of the scale, . For convenience, I will subsequently suppress 
the dependence on Q^. The asymmetry is given by the difference in probability 
of finding a gluon with its polarization parallel versus antiparallel 
to the proton's longitudinal polarization. The impolarizcd gluon distribution 
function is given by G{x) = G^{x) + G~(x). The gluon polarization is given 
by the ratio, AG{x)/G{x). Similar definitions of unpolarized parton distribu- 
tion functions and helicity asymmetry distributions exist for quarks (q) and 
antiquarks (g). 

The importance of AG is twofold. First, the proton's spin can be decom- 
posed as follows: 

5, = i = ^AS + AG + L« + Lf. (3) 

Here, Li^^^ gives the orbital angular momentum contributions of quarks (glu- 
ons) to the proton's spin. The contribution quarks make to the proton's spin 
(AE) is determined from quark helicity asymmetry functions, Aqi{x), summed 
over the quark flavors (n/) consistent with Q^, by: 

AS = ^Agi, where Aqi = / [q+{x) +qf{x) - q7{x) -qr{x)]dx. (4) 
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Hence, if AE is known, then a determination of AG establishes the contribution 
of partonic orbital motion to the proton's spin. 

The second reason for the importance of AG is that, due to the axial 
anomaly of QCD, the quantity AE cannot be determined independently of AG 
in polarized deep inelastic scattering (PDIS). This is graphically illustrated 
in Fig. 2, showing the results from the global anabasis of polarized deep- 
inelastic scattering (PDIS) made by the SMC gmupEI. The analysis of the 
limited information on scaling violations in PDIS 113 has produced only crude 
constraints on AG(a;) and its integral. 

The question then is, how can AG be determined? There are several 
possible methods being pursued: 

• Measurements of PDIS spanning a broad range of x and could be 
performed. A determination of AG(a::) would result from analysis of the 
scaling violations from this more extensive data set. This method would 
require a high-energy polarized ep collider. Possible plan&_£pr pursuing 
such a program at HERA were discussed at this workshop til. 

• Measurement of the leptoproduction of di-jet events is sensitive to AG(a::) 
through the photon-gluon fusion process. A variation of this method is 
to detect high-p^ charged hadron pairs, assumed to be the leading par- 
ticles of jets, in high-energy polarized-lepton/polarized-proton collisions. 
This method is being pursued by COMPASS 113, and intriguing data has 
already been obtained at much lower energy by HERMES t3. 
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• Measurements of di-jet production, high-pr particle or photon produc- 
tion in p + p collisions arc sensitive to AG{x). This is the primary focus 
of the RHIC-spin program. 

Sensitivity to AG{x) in polarized proton collisions arises for partonic col- 
lisions involving gluons. Di-jct production can be initiated either by qq, qg 
or gg collisions, with the latter two possibilities involving AG{x) linearly or 
quadratically, respectively. Just as for lepton-induced processes, detection of 
high-pT hadrons, assumed to be the particles of final-state jets, is also sensitive 
to AG(x). But, the relationship between the initial-state partonic kinemat- 
ics and the hadron's px is more complex because of the momentum sharing 
between the multiple hadrons within the jets. 

Direct photon production is, in principle, the cleanest probe of AG{x) 
for p -hp collisions, since in leading-order pQCD, observables are only linearly 
dependent on the gluon structure function for the dominant gluon Compton 
scattering process gq jq. For pp "fX. there is only a small physics 
background from qq 75. Furthermore, when the photon is detected in 
coincidence with the away-side jet, the initial-state partonic kinematics can be 
reconstructed (sec Sect. 5.3). Hence, direct photon production is a primary 
focus of the RHIC spin program. 

5 Plans for determining AG at STAR 

One of the most promising methods for determining AG is to study direct 
photon production in polarized proton collisions. Both PHENIX and STAR 
will study this process for the RHIC-spin program. Below, I give a detailed 
assessment, somewhat biased towards the performance of the STAR detector, 
of how accurately AG will be determined in polarized proton collisions at 
RHIC. 

5. 1 General features of direct photon production in pp collisions 

Within the framework of leading-order pcrturbative QCD, the dominant mech- 
anism for producing a single photon with large transverse momentum in a pp 
collision is gluon Compton scattering {gq jq) . The competing partonic sub- 
process is qq annihilation {qq jg). Given the predominance of gluons over 
antiquarks at small Bjorken x, gluon Compton scattering is ten times more 
likely to occur. Fig. 3 shows the pQCD predictions for the angular depen- 
dence of the cross section and longitudinal spin correlation coefficient for these 
two processes in the partonic center-of-mass reference frame. The polarized 
cross section (Act = cillct) is strongly peaked at scattering angles where the 



7 



1 







-1 



1 




_ 


1 







photon is emitted in the direction of the incident quark. This result is inde- 
pendent of the total energy in the partonic CM frame. The best determination 
of the gluon polarization will be made when the final-state photon is parallel 
to the initial-state quark; at other angles, there is reduced sensitivity to the 
gluon polarization. 

A second important criterion for opti- 
mizing the sensitivity to the gluon polariza- 
tion is to require large polarizatipns for the 
initial-state quark. Measurements L3 of polar- 
ized deep-inelastic scattering (DIS) have de- 
termined that the quark polarization increases 
linearly with log x, for momentum fractions 
greater than 0.1. At « 0.2, the quark po- 
larization is ~ 0.3, and it continues to rise as 
Xq increases. Furthermore, we have learned 
from unpolarized DIS, that the small x region 
{xg < 0.1) iSpflrecisely where gluons predomi- 
nately resideEil. Hence, the optimum determi- 
nation of the fraction of the proton's spin car- 
ried by gluons will be made by using large-x 
quarks {xq > 0.2) as an analyzer of the polar- 
ization of small- a; gluons. 

These two criteria suggest that the great- 
est sensitivity to AG(a;) results from asymmet- 
ric gluon Compton scattering (xq > Xg), with 
the photon detected in the same direction that 
the partonic CM moves in the collider refer- 
ence frame. As a consequence, both the pho- 
ton and the hadronic jet from gluon Compton 
scattering events should be detected at large 
pseudorapidity, limited only by the need for 

large pr in the collision, generally assumed to ative to the incident gluon. 
be proportional to the kinematic scale, Qi relevant for the structure functions. 
To limit the contributions from higher-twist processes, we will consider below 
direct photon production with pt,^ > 10 GeV/c. 




Figure 3: Leading-order perturba- 
tive QCD predictions for the par- 
tonic spin correlation coefficient 
and the relative differential cross 
section for photon production pro- 
cesses. For gluon Compton scat- 
tering, the partonic center of mass 
scattering angle (d*) is defined rel- 



5.2 Simulating polarized proton collisions 

Sophisticated event generators have been developed for high-energy (unpolar- 
ized) pp and pp collisions to aid in the optimization of experiments. To explore 
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questions about the expected perforraajice of the STAR detector we turn to a 
QCD event gener ator (PYTHIA 5.7)113, known to include many of the salient 
features of hard scattering processes. Hadronization of the recoiling quarks and 
gluons following a hard-scattering event is accounted for by the Lund string 
model, with parameters tuned to agree with fragmentation functions measured 
in e+e~ colliders. The multiple soft-gluon emission, thought to be responsible 
for the introduction of transverse momentum to the initiaLstate partons, is 
accounted for in PYTHIA through the parton-shower modelll3, rather than by 
explicit evaluation of higher-order QCD processes. 

Spin effects are included for each event by evaluating the leading-order 
pQCD expressions for the process-specific spin-correlation coefficient 113, a^™'^, 
using the Mandelstam variables {s,t,u) for the partonic hard scattering, as 
given by PYTHIA. For the polarization |Observables, the helicity asymmetry 
distributions from Gehrmann and Stirlingll3 are evaluated at the Bjorken x val- 
ues (a;i(2)) given by PYTHIA, after being evolvedllS to = Gehrmann 
and Stirling provide three different sets of distributions, referred to here as GS- 
A,B and C. The sets differ mostly for AG{x) (Fig. 2). Armed with these vari- 
ables, the proton spin correlation for a given pp scattering event is calculated 
as: 



^LL = , , . / r,7\ O-LL {S,t,u). (5) 



In this equation, Af / f is the polarization of the parton (either a gluon or 
a quark (antiquark) of a given flavor) within one of the interacting protons 
(either beam a or 6), assuming the beams are 100% polarized. The generated 
event may in principle arise from colliding protons with either equal (-I-) or 
opposite (-) helicities, but with different probabilities, proportional to 

^i± = [l±Pb,Pb,ALL\(Jeff- (6) 

In the simulations, the initial helicity state for a given event is chosen ran- 
domly in accordance with the above relative probabilities. In Eqn. 6 the beam 
polarizations, Pbn^)! ^re each taken to be equal to 0.7, as expected for the 
RHIC-spin program. Alternate selection between the two polarization states 
(representing for two-spin observables either equal or opposite helicity states 
for the colliding protons) is continued until a non-zero value is drawn, corre- 
sponding to the occurrence of a collision. 

A key assumption in this method of computing polarization observables, 
is that initial-state parton showers are snin independent. This assumption 
has been checked in simulations (SPHINX E3) that separately consider parton 
showers independently for each helicity projection- of the initial-state partons 
before the hard scattering. The good agreement E2l between SPHINX and the 
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Figure 4: Simulated values for All for inclusive photon production at RHIC energies. The 
top row shows the spin correlations for midrapidity photons that can be detected at both 
PHENIX and STAR. The bottom row shows expected values at the forward angles probed 
by the STAR endcap electromagnetic calorimeter (EEMC). Not evident in this figure, is that 
for a given xt, photons detected in the EEMC correspond to smaller-x gluons than those 
detected at midrapidity. 



method outlined above for the polarization observables suggests minimal spin- 
dependent effects from the initial-state parton showers. Further siipport for 
this conclusion comes from next-to-leading order QCD calculations EiJ. Higher- 
order processes are found to minimally influence the polarization observables 
computed in leading order. One precaution is that, to date, the necessary 
resummation, required to account for multiple soft gluon emission, has not yet 
been carried out for photon production processes. 

Given this methodology, polarization observables have been calculated for 
a variety of processes. The acceptance and finite resolution of the STAR barrel 
and endcap calorimeters (F, MC) for photons is accounted for in these simula- 
tions. Isolation criteria l3 are applied to the simulations, anticipating their 
need in real experiments to reduce backgrounds (Sect. 5.5). Some of the 
calculations presented below rely on hadronic jet reconstruction. For those 
calculations, the response of the STAR TPC and the EMC are approximately 
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accounted for. Standard jet reconstructfon algorithmsE3 are applied, based on 
the charged particle tracking from the STAR TPC and the electromagnetic 
energy detection from the EMC, and assuming that long-lived neutral hadrons 
{n,n,Kj^) are not detected. To avoid edge effects, the reconstructed jet is re- 
quired to be within a full cone radius (taken to be of size, R = 0.7 radians) 
from the edge of the STAR detector hmits. 

The resulting asymmetries for inclusive direct photon production, under 
the conditions described above, are shown in Fig. 4 as a function of the photon 
transverse momentum, scaled by ^/s. Several observations can be made: 

• For the GS-A,B gluon helicity asymmetry distributions, the All values 
are expected to be 'large' (> 0.05) for large xt, corresponding to values 
of Xg where the gluon polarization is large. The small All values for 
GS-C reflect the small gluon polarization at all Xg, arising because the 
peak of xAG{x) (Fig. 2) occurs at an Xg value where the unpolarized 
distribution is already quite large. 

• The qualitative features of direct photon production (discussed in Sect. 
5.1) are evident by comparing All for midrapidity and forward-angle 
photons. The latter, in general, results in larger magnitude asymmetries; 
and hence, a greater sensitivity to AG. 

• It is necessary to detect very hard photons {pr,-, > 30 GeV/c) to enable 
measurements All at the same xt for two values of ^/s. Higher 
energy photon detection is more difficult, but is important to accomplish. 
Sampling the same xt at different ^/s is an important check of the theory. 
Direct photon production cross sections have been shown to deviate from 
next-to- leading order QCD calculations in an xr-dependcnt manner O. 
One speculation for this is the need for larger kx smearing than is present 
in the theory. 

5.3 Determining the initial-state partonic kinematics 

Direct photon production, where only the photon is detected (inclusive detec- 
tion), provides only crude determination of the initial-state partonic kinemat- 
ics. It is generally assumed for photons detected at mid-rapidity (77 « 0), that 
the kinematic quantity xt — '2pT,-y/\/s is approximately equal to the Bjorken 
X values of the initial-state colliding partons. This assumption is only approxi- 
mately valid for the average values of Xg and Xq. In fact, due to the lack of any 
constraint on the direction of the recoiling jet, Xg and Xq cirG distributed along 
an approximately hyperbolic locus at fixed xt and rj-y. This locus is slightly 
broadened by kx smearing effects. 
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reconstructed J^^j,^,.^. reconstructed x^,,,^,. 

Figure 5: Correlation between reconstructed and simulated values for the quark and gluon 
momentum fractions for the p + p— > ■y + jet + X reaction. The simulations include effects of 
kx smearing. The reconstruction assumes coUinear collisions. Only events with an isolated 
photon in the STAR acceptance (— 1 < »7 < 2) having 10 < px < 20 GeV/c are accepted. 
The coincident jet must be within the fiducial volume of STAR (—0.3 < Tj^et < 1.3). A 
further restriction that the reconstructed Xq^^rk > 0.2 is also imposed. 



These kinematic ambiguities can be vastly reduced by detecting the away- 
side jet in coincidence with the photon. With the reasonable assumption of 
collinearity of the colliding partons in the initial state, it is easy to show that 
the initial-state momentum fractions can be determined for each event by mea- 
suring the energy {pT.-t) and direction of the photon (77^,), and only the direc- 
tion of the away-side hadronic jet {r]jet)- With this information, conservation 
of energy and momentum at the partonic level imply 

xi = ^(e-"^ + e-"^"') and X2 = ^(6+""' + e+^^^O- (7) 

Given the measured quark and gluon helicity-independent probability distri- 
butions, the reasonable assumption is made that the initial-state quark had 
momentum fraction a;^*^™" = max[a;i,a;2] and the gluon, Xg'^'^°"' = min[2;i, 0:2]. 
Similarly, it is straightforward to express the partonic CM scattering angle in 
terms of rjj and rjjet- A key assumption, that is examined more carefully be- 
low, is that the initial-state partons are coUinear. A valid question is, to what 
extent do transverse momentum components in the initial state (fcr) interfere 
with the event-by-event determination of the kinematics? 

To address this question, the reconstruction algorithm was applied to 
PYTHIA simulations of direct photon processes. The initial-state partons had 
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transverse momentum components, as generated by the parton shower modeL 
Events having at least one jet and a coincident photon within-thc STAR fidu- 
cial volume were reconstructed. The UA2 isolation conditioner was applied to 
the photon. A condition that the reconstructed Xq was greater than 0.2 was 
also imposed. As is evident in Fig. 5, most of the events have the initial-state 
momentum fractions properly reconstructed. The majority of the kinematic 
reconstruction errors occur when Xq < Xg. 

5.4 Direct reconstruction of AG{x) 

Armed with the knowledge that the initial-state partonic kinematics can be ac- 
curately reconstructed, it is possible to consider directly reconstructing AG{x) 
from the measured longitudinal spin correlation. It is possible that such a 
direct reconstruction will provide the best framework for deducing AG from 
photon production measurements. As is already evident from Eqn. 5, if only 
quark-gluon Compton scattering contributes to the photon yield, then 



All = ^GM2^ A\{xq^ Q2)a^r'*°"(^*)- (8) 



This equation results from Eqn. 5, when averaging over an event ensemble, 
since the u and d quark contributions to gluon Compton scattering are weighted 
by their squared electric charge and the probability to find them inside the 
proton. This is identical to the weighting that enters in polarized deep-inelastic 
scattering (PDIS).|Xhe quantity, A^, in Eqn. 8 is then precisely the asymmetry 
measured in PDISEJ. 

Since the kinematics can be determined for individual events, it is possi- 
ble to invert Eqn. 8 to express the reconstructed gluon helicity asymmetry 
distribution in terms of known quantities: the unpolarized gluon PDF, the 
quark polarization and the pQCD expression for the gluon Compton scat- 
tering partonic spin correlation; and the measured spin-dependent scattering 
yields (A^++(+^)): 

. ^ [N++{Xg) - N+^{Xg)] 

(9) 

This inversion, as applied to simulations oi All for direct photon processes is 
shown in Fig. 6. The key assumptions implicit to this direct reconstruction, 
and their influence on /S.Grecon{x), are: 

• Only the gluon Compton process contributes to the 7 -/-jet yield; the other 
processes that contribute to the simulated yields are assumed to be absent 
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Figure 6: Results for the direct reconstruction of AG(x,Q^), as applied to the simulated 
values of Ai^j^ for direct photon processes. Three different inputs for AG(x) are shown. 



in the reconstruction. The small contribution from qq — s- 717 is negative, 
and is proportional to the product of the quark and antiquark polar- 
izations at the asymmetric Bjorkcn x values probed in the experiment. 
This contribution results in AGrecon{x) underestimating the input value 
of AG(x). 

• The initial-state partonic kinematics are assumed to he perfectly recon- 
structed. For some events, especially the small number with Xg > Xq, 
there are reconstruction errors. Typically, All for the Xg value recon- 
structed for these events is smaller than for those events that don't have 
kinematic reconstruction errors. The end result is that AGreconix) un- 
derestimates the input value of AG{x). 

• The partonic collisions are assumed to be collinear, so kx smearing effects 
are assumed to he ahsent. These effects introduce non-zero transverse 
momentum in the initial state, contrary to the assumption of collinear 
collisions. The end result is to make correlated errors in the reconstruc- 
tion of Xq and Xg. 

Corrections to AGrecon{x) for these three effects can be made based on 
simulations. The first correction requires no knowledge of AG{x), whereas the 
latter two would need to be made in an iterative fashion, since knowledge of 
AG{x) is required. 
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In addition to these three effects, present in the results shown in Fig. 6, 
other effects can influence AGreconix) that will be determined from p + p 
J + jet + X data. In particular, the 77 decay of a high-energy 7r°(77) mesop can 
mimic the detector response of a single high-energy photon. Simulations^ have 
shown that this background results in a smaller magnitude All compared to 
direct photon production, and if not properly corrected, would result in an un- 
derestimate of the true AG{x) via the direct reconstruction method. Another 
problem facing the RHIC-spin experiments, is the possibility that some frac- 
tion of the photons produced in p -I- p collisions arise from the 'fragmentation' 
of recoiling final state partons. In particular, high-pT photons can be produced 
in the hard bremsstrahlung of a charged parton, produced by qg or gg scatter- 
ing processes that havej-significantly larger cross section than the direct photon 
processes. SimulationsEa have shown that these other processes will dilute the 
direct photon All by a small amount. Again, they tend to underestimate the 
true AG(x) via the direct reconstruction method. 



5.5 Backgrounds and their suppression 

To employ direct photon production as a means of determining the polarization 
of the proton's glue, a detector must be capable of selecting the very small 
fraction of events that have a single high energy photon. The most pernicious 
physics background arises from di-jet events (by far, the largest fraction of the 
total reaction cross section), where one of the jets has an energetic neutral 
meson (tt^ or 77°, collectively referred to as M") that can decay into a pair of 
photons. Kinematically, the most probable opening angle between the photons 
produced by M° decay is 0™'' = 2sin-^{mMoc^ / Emo)- For a 30 GeV 7r° this 
angle is 9 mr. Such a small opening angle makes it very difficult for any 
detector to distinguish between di-photons from AI^ decay and a single direct 
photon. 

The relative probability for producing a neutral meson versus a direct 
photon is shown in Fig. 7a. Due to the larger number of processes that can 
produce jets versus single photons, and the fact that the strong interaction 
coupling constant (as) is involved rather than the electromagnetic coupling 
constant (a), Af" production is nearly an order of magnitude larger than 7 
production. This result from PYTHIA has been shown tJ to be roughly in 
accord with existing measurements E^. 

The influence of the background on the direct photon measurements 
at RHIC is different for the PHENIX and STAR detectors. The former has a 
smaller solid angle, but substantially finer granularity, EMC than the latter. 
This feature allows PHENIX to better reconstruct M° from their daughter 
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photons. The uniform distribution of the decay photon energies from de- 
cay ultimately sets a limit to the range of decay phase space where this back- 
ground suppression technique is effective. The STAR detector will distinguish 
background from signal using two primary techniques, described below. 
The first of the background sup- 



pression methods for STAR relies 
on a 'shower maximum detector' 
(SMD^. This is a fine granularity 
detector placed within the depth of 
the calorimeter to measure the trans- 
verse profile of the electromagnetic 
(EM) showers produced by the in- 
cident particles. In general, an in- 
cident M° results in two separated 
clumps of energy deposition in the 
SMD, corresponding to the spatial 
separation of the EM showers pro- 
duced by the closely spaced daugh- 
ter photons. An incident direct pho- 
ton typically produces only a single 
clump- of energy deposition in the 
SMD I. This distinction is mud- 
died by the complexity of EM show- 
ers, resulting in false dumpiness in 
the SMD response for some fraction 
of the direct-photon-induced events. 
Despite this, the analysis of the SMD 
information will provide significant 
suppression of the M° background. 
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Figure 7: Ratio of cross sections for n^{ri) and 
direct photon production under increasingly 
stringent conditions: a) including only a 
cut; b) cuts on kinematic quantities and an 
isolation cut; c) cut on information from SMD. 

The second component of the background suppression arsenal of STAR is 
related to the nearly 47r coverage the detector provides. In general, high energy 
neutral mesons are leading particles of jets, and hence, are accompanied by 
additional hadrons. An effective discrimination between M'^ and direct photon 
events is provided by examining the event topology. The former class of events 
generally have additional hadrons within a relatively narrow cone around the 

whereas direct photon events have little accompanying energy within the 
isolation cone. The M" events can be suppressed by using 'isolation cuts' 
similar to those used in other collider detectors reporting direct photon cross 
sections. 

The influence of these background reduction cuts on the background:signal 
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ratio is illustrated in Figs. 7b and 7c. The end result is that the signal is ex- 
pected to be approximately 2x larger than the background. The remaining 
background contributions can be eliminated by taking the difference between 
distributions enriched in direct photons and those enriched in M'^. This sub- 
traction process will increase the statistical errors beyond those shown in Fig. 
6 by a factor varying between 1.5 and 2u. The systematic error associated 
with these corrections can be estimated from in situ calibrations of the perfor- 
mance of the SMD, made possible by producing an energetic tt^ sample from 
reconstructed events. Other backgrounds will also be present in both the 
STAR and PHENIX experiments, but their effects on the determination of AG 
should be less important. Efforts are underway to understand the limitations 
on extracting AG from direct photon measurements. 



5.6 Extrapolation Errors 

Unlike the goals_Df most experiments sensitive to gluon polarization that are 
either underwayEj or are on the horizonlij, the RHIC-spin program aims to de- 
termine the integral contribution gluons make to the proton's spin, as defined 
in Eqn. 1. To carry out this integral, the gluon helicity asymmetry distri- 
bution must be determined over a sufficiently broad range of Xg to minimize 
the influence of extrapolation errors in carrying out the integral over all x. 
To illustrate how these errors influence the determination of AG, a standard 
parameterization of the Xg dependence, 

xAG{x) = r^Ax^il ~ xf[l + px^/^ + jx], 
V a + b+lj T{a + b + l) ^ r(a + 6+l) ^ ' 



so that rj = AG{x)dx, 
Jo 



is fit to simulated data. Similar to the method used to analyze existing data 
on scaling violations in polarized deepJLnelastic scattering E3, b and 7 are held 
fixed at values obtained by evolving 13 the AG(x) input to the simulation 
to the values relevant at RHIC. The fixed parameters used in the fits 
are consistent with positivity constraints (|AG(a;)| < G{x)). The parameters 
?7, a, p are then adjusted to provide the best fit to AGrecon{x), described in Sect. 
5.4. The results for the simulations of the y/s—200 and 500 GeV samples are 
combined, after making additive corrections for qq annihilation contributions 
to the simulated asymmetries. No corrections have been made for kinematic 
reconstruction errors or kr smearing. As well, the values of AGrecon{x) have 
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Figure 8: Fits to the reconstructed AG(a:), after correcting for qq annihilation. A standard 
parameterization (Eqn. 10) is used to fit the data, with the b and 7 parameters fixed. 
Full data sets for \/S=200(500) GeV are assumed, corresponding to integrated luminosity of 
320(800) pb-i. 



not been evolved to a common . These corrections are important, but have 
not been made because they require knowledge of AG'(x). The fits are shown 
in Fig. 8. It is found that: 

• an accurate determination of AG will require both ^/s—200 and 500 GeV 
data samples to get to sufficiently small Xg. Due to strong correlations 
between r} and a, 6ri grows rapidly as the low- a; points are successively 
eliminated. The error in the integral AG (^77), wijtb both samples in- 
cluded in the fit, is 0.24 when Gehrmann-Stirling t2l set A is input to 
the simulations. The value for 5rj for set B is 0.62, and the relative er- 
ror in ?7J£ comparable for set C. For set A, other systematic errors are 
expecteoH to increase the error in the integral AG to 0.5, still sufficient 
to accurately establish AE from polarized deep inelastic scattering. 

• the large-x behavior of AG(a;) must be constrained to determine AG at 
RHIC. Different values for the fixed parameters {b, 7) result in different 
values for 77. 

• fitting /S.Grecon{x), including only the corrections for qq annihilation, 
yields a value for the fitted rj that is too small compared to the input AG. 
Evolving all of the AGreconix) points to a common and correcting 
for the kinematic reconstruction errors results in a fitted 77 that is closer 
to the input AG, but is still too small. Both of these corrections require 
knowledge of AG(x), and hence will require an iterative procedure for 
their application. The largest remaining error comes from kx smearing. 
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Repeating the analysis with simulations that don't include initial-state 
parton showers results in a fitted 77 in agreement with the input AG. 

The end result_is, that after accounting for the most significant sources 
of systematic error we expect that the fraction of the proton's spin carried 
by gluons can be determined to an accuracy of approximately 0.5, primarily 
based on the STAR measurements oi p + p 7 + jet + X. Data samples at 
both -^5=200 and 500 GeV are crucial so that the accuracy is not limited by 
extrapolation errors. The analysis of AGrecon{x) presented here is intended 
to illustrate the sensitivity of the STAR measurements to the integral AG. 
Clearly, the best determination of AG will result from a global analysis of all 
relevant data. 



6 Summary 

The RHIC-spin program promises to provide exciting data to test perturba- 
tive QCD, and to provide new insights into the non-perturbative structure of 
the proton. It is likely that one of the most interesting results will be the 
determination of the fraction of the proton's spin carried by gluons. I ex- 
pect that inclusion of AGrecon{x), based on the measurement of All for the 
p + p 'j + jet + X reaction at -^5=200 and 500 GeV, in a global analysis will 
provide the best determination of the gluon contribution to the proton's spin. 
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